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Abstract 

A computational hemodynamics approach is presented to simulate the blood flow through the human circulatory 
system under altered gravity conditions. Numerical techniques relevant to hemodynamics issues are introduced to 
non-Newtonian modeling for flow characteristics governed by red blood cells, distensible wall motion due to the heart 
pulse, and capillary bed modeling for outflow boundary conditions. Gravitational body force terms are added to the 
Navier-Stokes equations to study the effects of gravity on internal flows. Six-type gravity benchmark problems are 
originally presented to provide the fundamental understanding of gravitational effects on the human circulatory 
system. For code validation, computed results are compared with steady and unsteady experimental data for non- 
Newtonian flows in a carotid bifurcation model and a curved circular tube, respectively. This computational approach 
is then applied to the blood circulation in the human brain as a target problem. A three-dimensional, idealized Circle 
of Willis configuration is developed with minor arteries truncated based on anatomical data. Demonstrated is not only 
the mechanism of the collateral circulation but also the effects of gravity on the distensible wall motion and resultant 
flow patterns. 


Introduction 

Currently, NASA has been challenging long-duration 
space missions for the International Space Station (ISS) 
and flight to Mars. During the elongated space mission, 
astronauts have to adapt themselves to the altered cir- 
cumstance of microgravity. Blood circulation as well as 
body fluids distribution undergoes significant adaptation 
during and after space flight. Much study on physiologi- 
cal changes under weightlessness has been performed 
since the early days of the space program [1]. In particu- 
lar, cardiovascular research in conjunction with the 
Space Shuttle program has included diverse physiologi- 
cal functions affected by the nervous system such as 
heart rate, blood pressure, hormone release, and respira- 
tion. The altered cardiac output due to adaptation during 
flight and deconditioning after the flight will impact 
blood circulation in the human body. Especially, this 
altered blood supply in the brain and consequent oxygen 
supply to certain parts of the brain will make non- 
negligible impact on long-duration flight. Also gravity 
over 8G causes unconsciousness, so-called, blackout 
whereas gravity below -3G makes the retina engorged 
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with blood, red-out. To assess the impact of changing 
gravitational forces on human space and earth flight, it 
will be essential to quantify the flow characteristics in the 
brain under varying gravity conditions. Analysis of blood 
circulation in brain as well as other parts of human body, 
such as heart, requires the capability to analyze flow in 
large arteries and capillaries. In addition to the altered 
gravitational forces, specific shapes and connections of 
arteries in the brain vary in the human population [2, 3]. 
Considering the geometric variations, pulsatile unsteadi- 
ness, and moving walls, computational approach in ana- 
lyzing altered blood circulation will thus offer an eco- 
nomical alternative to experiments. 

Major technical challenges of the present study relevant 
to hemodynamics issues can be summarized as i) non- 
Newtonian modeling of the human blood, ii) distensible 
wall motion due to heart pulse, iii) capillary bed modeling 
for outflow boundary conditions, iv) anatomically com- 
plex geometry of vascular networks, and v) gravitational 
effects. Regarding the non-Newtonian modeling, blood 
flow characteristics are dominated by the behavior of red 
blood cells since red blood cells constitute the human 
blood about 40 to 45 percent in volume. Red blood cells 
aggregation at low shear rates makes apparent blood vis- 
cosity rapidly increased. Conversely, the viscosity de- 
creases slowly at higher shear rates with red blood cell 
deformability. Three different experiments [4-6] in vitro 
show this so-called shear-thinning phenomenon in Fig. 1. 
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Various non-Newtonian models have been developed by 
curve fitting to the experimental data [7]. They assume in 
common that apparent blood viscosity is a function of 
shear rate alone. In the present study, two different mod- 
els are assessed through steady and unsteady non- 
Newtonian flow simulations. More details about both 
models are found in Computational Approach. 

Relative diameter change due to the heart pulse is 
found to be up to about 10 percent in common carotid 
arteries of young people [8]. Normal arteries have elastic 
lamina to transfer the blood most efficiently whereas 
veins have no elastic lamina but valves to prevent re- 
verse flows. Both arteries and veins in vivo show de- 
formable wall motions according to the supine or erect 
posture under certain gravity. On Earth, gravity pulls 
down the blood to the feet. In standing posture, the blood 
pressure in the feet can be high up to 200 mmHg (milli- 
meters of mercury) whereas it is about 60 mmHg in the 
brain. In space, blood pressure equalizes and becomes 
about 100 mmHg throughout the human body. Conse- 
quently, the arterial wall distensibility due to gravity 
causes fluid redistribution throughout the whole body. In 
the present study, a standard first-order approximation [9, 
10] to the complex behavior of the arterial wall is 
adopted by assuming that blood vessels have circular, 
thin and elastic walls. 

To accomplish the simulation at a computationally 
manageable level, minor arteries such as arterioles, 
venules, and capillaries need to be truncated. At the trun- 
cated positions where flow information is hardly obtain- 
able in vivo, equivalent outflow boundary conditions are 
required. In this study, a novel approach utilizing the 
analogy between arterial networks and electric circuits 
based on Poiseuille’s formula is introduced as an alterna- 
tive to impose adequate outflow conditions. The alterna- 
tive approaches using the analogy of electric circuits to 
arterial networks have been presented to supply the 
boundary conditions for three-dimensional computations 
[ 11 , 12 ]. 

For anatomically accurate three-dimensional recon- 
struction of complex arterial configuration, image seg- 
mentation techniques have been presented using mag- 
netic resonance angiogram (MRA) [12], magnetic reso- 
nance imaging (MRI) [13], and computed tomography 
(CT) [14], There have been some remarkable studies 
about computational simulation of the brain circulation 
[12,15,16]. They assumed the flow Newtonian in large 
arteries with solid walls. To the authors’ knowledge, the 
present study is the first one to demonstrate the blood 
circulation in the human brain under altered gravity con- 
ditions, by assuming non-Newtonian flow within disten- 
sible arterial walls. 


Computational Approach 
Governing Equations 

The flow through the heart and blood vessels is un- 
steady, viscous, and incompressible. The blood flow 
through heart and large vessels has been assumed Newto- 
nian and incompressible because blood is a liquid with a 
viscosity about 3.5 times that of water and acts as a simple 
fluid with suspended particles such as red blood cells [12- 
18]. For instance, the Newtonian flows through artificial 
heart devices have been simulated using the INS3D code 
[17, 18]. In the present study, the flows of interest are as- 
sumed non-Newtonian to account for the shear thinning 
behavior of the human blood. Governing equations are the 
three-dimensional, unsteady, incompressible Navier- 
Stokes equations coupled with non-Newtonian models, 
and can be written in tensor notation form as: 
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The shear stress tensor, r t j is defined as: 
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where M and rj are molecular and apparent viscosity, 
respectively. Sy is the mean strain-rate tensor and y is 

the shear rate defined as a function of the second invariant 
of S t j in three-dimensional problems: 
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Also the body force terms of gravity, g t are added to 
Eq.(2) to take gravitational effects into account. 

Non-Newtonian models 

Two different types of non-Newtonian models are im- 
plemented into the INS3D code and assessed through 
steady and unsteady computations. One is a Carreu- 
Yasuda model and the other is an extended Casson model. 
These models have been used in some computational 
hemodynamic simulations [ 1 9-21]. 

A Carreau-Yasuda model in Eq. (6) describes the shear 
thinning behavior of blood flows with asymptotic appar- 
ent viscosities at zero and infinite shear rates: 
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- n — =[l+(#) a f^ (6) 

Vo ~V°° 

where X and a define the width of the transition region 
from Newtonian to power-law region, and n - 1 is power- 
law slope as shown in Fig. 1. The constitutive parameters 
for the human blood are given by: 

r] „ = 0.00348Pa • s, rj Q = 0.1518Pa • 5 , 

' /0 (7) 

X = 40.05, a - 2.0, n = 0.356 

For the extended Casson model, 

7 = (-[c I (//r)+C 2 (/ffX/7f.V max ) (8) 

Y 


where Ht is the hematocrit, Cj and C 2 are coefficients 
determined for Ht = 40 percent as 

C, =Q.2(dyn/cm 2 ) in and C 2 = 0.18(<fyn- j/cm 2 ) 1 ' 2 
based on the experimental data [5]. To avoid extreme 
values at lower shear rates, the apparent viscosity is con- 
fined within /7 max =0.1518Pa*5. 

Deformable wall model 

A first-order approximate model for arterial wall dis- 
tensibility is adopted based on the assumption of a thin- 
walled and linearly elastic channel. Although not a com- 
plete model of the arterial wall, this approach was found 
to be a standard first-order approximation to the complex 
behavior of the arterial wall [9, 10]. Tethering to sur- 
rounding structures is thought to reduce the arterial 
longitudinal motion to a minimum [22]. In the present 
study, blood vessels are assumed to have circular, thin 
and elastic walls with negligible longitudinal motion. 
The increments of blood vessel radius are given by: 

-y- = — -Pwh =D w {p w -p w ) t (9) 


where p w and p w are transmural and reference pres- 
sure at wall, and the arterial wall distensibility D w is 
defined as: 
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where h is the wall thickness, E and v p are the elastic 

modulus and Poisson’s ratio of the arterial wall, respec- 
tively. The arterial wall distensibility is given from 
9.45 ~ 26.4 xl0~ 7 cm 2 I dyne based on the structural 


wall properties [10, 23]. For deformable walls, non-slip 
boundary conditions are replaced by moving wall bound- 
ary conditions using the Eq. (9) at each physical time step. 

Capillary bed modeling 

To make the problem computationally manageable, mi- 
nor arteries such as arterioles, venules and capillaries need 
to be truncated. Equivalent outflow boundary conditions 
should be imposed at the truncated positions. There is an 
analogy between arterial networks and electric circuits as 
shown in Fig. 2. A truncated artery is assumed to divide 
into N branches of same size, for instance, N equals two 
for bifurcation and three for trifurcation. Under this as- 
sumption, the outflow boundary conditions, especially for 
pressure, can be approximately determined by utilizing 
this analogy and Poiseuille’s theorem: 

e = TT Ap=>e = ‘^'’ ^ + i=-^-Qk (ID 
SuL R H 


where Q is mass flow rate, p is pressure, and R is flow 
resistance. Also r and L are arterial radius and length, re- 
spectively. At the ^bifurcation (or trifurcation), its mass 
flow rate is N times of flow rate through the (£-1)* 
branches. Assuming that the flow resistance ratio /is con- 
stant, the pressure drop at each level can be expressed in a 
geometrical series form: 
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To avoid an unrealistic pressure difference between a 
truncated position and the capillary bed, the flow resis- 
tance ratio / should be determined less than one. For the 
problems with multiple outflow boundaries, it is prefera- 
bly given by 0.5 ~ 0.95. 

MPI-OpenMP parallel processing 

An MPI-OpenMP hybrid version of INS3D code [24] 
was extended with implementation of those numerical 
approaches, and validated through the comparison with 
experimental data for both steady and unsteady non- 
Newtonian cases. Computational experiments have been 
conducted on the INA Cluster (SGI workstation cluster) at 
NASA Ames Research Center. 
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Results and Discussion 
Steady non-Newtonian flow 

The steady non-Newtonian flow in a carotid bifurca- 
tion illustrated in Fig. 3 was simulated as the first valida- 
tion problem. Anatomically, common carotid arteries 
(CCA) in pairs arise from the ascending aortas of the 
heart and connect to the brain through the neck. CCA are 
divided into internal (ICA) and external carotid arteries 
(ECA), respectively. Figure 4 shows a chimera overset 
grid with 8 blocks. Computations are conducted with the 
same flow conditions of the experiment [19]. In this 
experiment, a blood analogy fluid (KSCN-X; KSCN- 
Xanthan gum solution) was used to mimic the shear 
thinning properties of blood. 

The constitutive parameters of the Carreau-Yasuda 
model in Eq. (7) are modified based on the experimental 
data as following: 

rj „ = 0.0022 Pa • s, rj 0 = Q.022Pa * r, 

(13) 

A = 0. 1 Is, a = 0.644, n = 0.392 

Likewise, a parameter of the extended Casson model is 
modified as C x = 0A(dyn/ cm 2 ) 112 for this blood anal- 
ogy fluid. 

The Reynolds number based on the CCA diameter is 
270, and flow division ratio of ECA over CCA is 0.45. 
The inflow through the CCA is assumed fully developed. 
The Figure 5 (a) and (b) show the comparison of axial 
velocity profiles in the ICA using two different non- 
Newtonian models. The results of the INS3D with chi- 
mera overset grids are compared with those of IFANS3D 
[25] with multiblock grids, and also compared to the 
experimental data. Axial velocity profiles are skewed 
towards the flow divider and the adverse pressure gradi- 
ent with the increase of the sectional area results in a 
secondary flow pattern in the sinus. The flow reacceler- 
ates after passing through the maximum diameter region. 
Overall both the computed results from two different 
models show an agreement with the experiment. Thus, 
the Carreau-Yasuda model was mainly used for the 
computations of non-Newtonian flows. 

Unsteady pulsatile non-Newtonian flow 

For an unsteady code validation, the pulsatile non- 
Newtonian flow in a 90 deg circular tube was simulated 
with the same flow conditions of the experiment [20]. 
The same blood analogy fluid with the steady case was 
used and the Reynolds number based on the tube diame- 
ter and diastole velocity is 300. The radii of tube and 
curvature are 4 and 24 millimeters, respectively. For pul- 
satile inflow boundary conditions, the experimental 
waveform was successfully regenerated using twelve 
harmonics based on the Fourier theorem as shown in Fig. 


6. Inflows are assumed fully developed and variant with 
time using the regenerated waveform. Computed results 
are compared with the experimental data at three different 
phases: end diastole, peak systole, and begin diastole. 
Figure 7 shows the comparison of axial velocity profiles 
on the plane of symmetry and the plane perpendicular to 
the symmetry plane. The computed results are coincident with 
the experimental data at all three phases. 

Gravity Benchmark Problems 

Six-type Gravity Benchmark Problems (GBP) are 
originally presented to provide the fundamental under- 
standing of gravitational effects on the human circulatory 
system. As shown in Fig.8, GBP-0 stands for microgravity 
whereas GBP-1 represents the human body in supine pos- 
ture under normal gravity on Earth. Among six problems, 
let us discuss GBP-2 and 3 cases representing the upper 
and lower part of the human body in standing posture, 
respectively. Test model is a U-tube with a diameter of 8 
millimeters and a curvature ratio of 4. The fluid is as- 
sumed Newtonian with a viscosity of 3.5cPoise and the 
Reynolds number based on the tube diameter is 169. Fig- 
ure 9 shows a comparison of axial velocity profiles be- 
tween solid- and deformable-wall cases at 90 deg station. 
A first-order approximation model in Eq. (9) is adopted 
for wall moving boundary conditions, and the wall disten- 
sibility factor is given about five times larger than that of 
normal arteries to bring gravitational effects into relief. 
Even though not all shown here, the gravitational effect on 
the internal flow within rigid wall is found to be minor. 


Table 1. Comparison of required power for GBP-2 and 3. 



Wall 

Ap 

(Pa) 

Q 

(cm 3 / s ) 

Ape 

(MW) 

GBP-2 

Rigid 

21.4151 

3.5186 

75.3512 

Deformable 

24.8279 

3.5186 

87.3595 

GBP-3 

Rigid 

21.0905 

3.5186 

74.2090 

Deformable 

19.1109 

3.5186 

67.2438 


For GBP-2W (GBP-2 with deformable wall), however, 
tube diameter is contracted by 5.1 %D because of the de- 
creased pressure due to gravity. Conversely for GBP-3W, 
the diameter is expanded by 5.3 %D because of the in- 
creased pressure. As flow resistance is, according to 
Poiseuille’s formula, inversely proportional to the fourth 
order of tube radius, GBP-2W requires about 16% more 
power whereas GBP-3 W requires about 10% less power 
as shown in Table 1. 

GBP provide a preliminary knowledge that gravity has 
considerable effects on deformable wall motion and con- 
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sequent flow patterns in the human circulatory system. 

Blood circulation in human brain 

There are two arteries in pairs that supply blood to the 
brain. One pair is common carotid artery and the other is 
vertebral artery. They are distally connected to the im- 
portant part of the brain, so-called, the Circle of Wills 
(CoW). The circle is sitting on the base of the brain and 
its main function is to distribute the blood evenly 
throughout whole of the brain. Regardless of race or age, 
its configuration is different among population. Actually 
it was found that only about 52 percent of the population 
have an ideal configuration [15, 16]. Thus, as a funda- 
mental approach, an idealized CoW geometry was de- 
signed based on anatomical data with minor arteries 
truncated, not only to demonstrate the gravitational ef- 
fects on hemodynamics in the brain, but also to simulate 
the mechanism of collateral circulation in the brain. Fig- 
ure 10 shows a chimera overset grid with 10 blocks for 
this ideally balanced configuration. Total number of grid 
mesh is about 300,000. 

Effects of altered gravity 
To demonstrate the effects of altered gravity, the com- 
puted result under the normal gravity of 1G was simu- 
lated and compared with those under microgravity or 
zero gravity. Figure 11 shows a non-Newtonian flow 
simulation through the balanced CoW configuration. 
Because of symmetry, blood flow through the anterior 
communicating artery (ACoA) is zero. Also the mass 
flux through the posterior communicating arteries 
(PCoA) is very low. For distensible wall motions due to 
the altered gravity, the wall distensibility factor, D w in 

Eqs.(9) and (10), was given by 26.4xl0“ 7 cm 2 J dyne . 

To avoid the unrealistic connectivity at the arterial con- 
junctions, the distensibility factor was set to be zero at 
each conjunction by assuming that wall distensibility at 
the conjunction is negligible compared with that of distal 
arteries. Distensibility factors at the transition region 
between the conjunction and arterial branches were in- 
terpolated using a smoothing function. The decreased 
blood pressure due to gravity results in arterial contrac- 
tion up to about 12 %D as zoomed in Fig. 12, which is 
quite similar to GBP-2W case among six-type GBR 

Collateral circulation 

Even though one of main arteries in the brain is 
stenosed or missing, the distal smaller arteries can re- 
ceive blood from the other arteries through the CoW. To 
simulate this interesting mechanism of “collateral circu- 
lation” in the brain, the left ICA is presumed occluded 
66.7 percent in flow rate. In other words, only one third 
of the blood is delivered through this artery to the CoW 


as indicated in Fig. 13. Unlike the balanced configuration 
case, the mass flux through the PCoA is most increased to 
compensate about 20 percent for the deficiency. Likewise, 
blood flow through ACoA is not zero but considerable. 
Instead, the mass flux through the proximal part of the left 
anterior cerebral artery (ACA) is very low in order to dis- 
tribute blood to the left middle cerebral artery (MCA) and 
left ACA simultaneously. This figure successfully demon- 
strates the mechanism of collateral circulation in the brain. 

* Further computational results and discussion will be 
presented in the final draft . 
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Shear Rate (1/s) 

Fig. 1 Logarithmic relation between blood viscosity and 
shear rate. 
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Fig. 8 Six-type gravity benchmark problems (GBP). 



(a) GBP-2 


(b) GBP-3 


Fig. 9 Comparison of axial velocity profiles between rigid- 
and distensible-wall cases. 


Fig. 7 Comparison of axial velocity profiles between compu- 
tation and experiment at three different phases. 
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(a) Microgravity (zero G) 


ACA: Anterior Cerebral Artery, ACoA: Anterior Communicating Artery, 
ICA. Internal Carotid Artery, MCA: Middle Cerebral Artery, 

PCA: Posterior Cerebral Artery, PCoA: Posterior Communicating Artery. 

Fig. 10 Design of an idealized CoW configuration. 


Fig. 11 Non-Newtonian flow simulation through the CoW 
under 1G 


(b) Normal gravity (1G) 

Fig. 12 Effects of gravity on arterial wall deformation and 
flow patterns of CoW. (zoomed view) 


Fig. 13 Collateral circulation through the CoW with left ICA 
66.7% stenosed. 
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